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Abstract 

We discuss a method to transform the covariant Fokker action into an implicit two-degree- 
of-freedom Hamiltonian for the electromagnetic two-body problem with arbitrary masses. This 
dynamical system appeared 100 years ago and it was popularized in the 1940's by the still in- 
complete Wheeler and Feynman program to quantize it as a means to overcome the divergencies 
of perturbative QED. Our finite-dimensional implicit Hamiltonian is closed and involves no series 
expansions. The Hamiltonian formalism is then used to motivate an EBK quantization based on 
the classical trajectories with a non-perturbative formula that predicts energies free of infinities. 
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I. INTRODUCTION 



The usual Hamiltonian description of two-body dynamics is surprisingly restrictive within 
relativity physics: If Lorentz transformations are to be represented by canonical transfor- 
mations, only non-interacting two-body motion can be described. This is the content of the 
no-interaction theorem of 1964, which later in 1984 was proved for the usual Lagrangian de- 
scription as well [1]. A covariant version of Hamiltonian dynamics, constraint dynamics, was 
invented to overcome this group-theoretical obstacle, but this too has limited applicability 
P). In this work we start from a relativistic physical theory: the time-symmetric relativistic 
actio n-at-a-distance electrodynamics, for which no constraint description exists at present. 
We discuss the passage from the covariant Fokker Lagrangian to an implicit Hamiltonian 
formalism defined in closed form without power expansions. We derive a two-degree of free- 
dom Hamiltonian for the electromagnetic two-body problem |3| with arbitrary masses and 
with either repulsive or attractive interaction. In our description a Lorentz transformation 
is represented by a canonical transformation followed by a rescaling of the evolution param- 
eter. We use this covariant Hamiltonian formalism to motivate an EBK quantization of the 
electromagnetic two-body problem that uses the classical trajectories to approximate the 
quantum energies with a formula that is non-perturbative and free of infinite quantities. 

In 1903, Schwarzchild proposed a relativistic interaction between charges that was time 
reversible precisely because it involved retarded and advanced interactions symmetrically 
[J]. The same model reappeared in the 1920s in the work of Tetrode and Fokker and it 
finally became electromagnetic theory after Wheeler and Feynman showed that this direct- 
interaction theory can describe all the classical electromagnetic phenomena (i.e. the classical 
laws of Coulomb, Faraday, Ampere, and Biot-Savart) jy, In particular, Wheeler and 
Feynman showed in 1945 that in the limit where the electron interacts with a completely 
absorbing universe, the response of this universe to the electron's field is equivalent to 
the local Lorentz-Dirac self-interaction theory [8] without the need of mass renormalization 



. Complete absorption is added to the theory as an approximation to uncouple away 
from the detailed neutral-delay dynamics of the other charges of the universe. For other 
approximations see Ref. [9j. The Wheeler and Feynman program to quantize the action- 
at-distance electrodynamics and overcome the infinities of QED is still not implemented 
because of the lack of a Hamiltonian description, which was a motivation for the present 
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work. 

The relativistic action-at-a-distance electrodynamics is better regarded as a many-body 
electromagnetic theory, once it is based on a parametrization-independent action involving 
two-body interactions only, without the mediation by fields jfj, [?J . The isolated electromag- 
netic two-body problem, away from the other charges of the universe, is a time-reversible 
dynamical system defined by the action 

Sf = — J m\ds\ — J m 2 <is2 — eie 2 / / 5(\ \xi — x 2 \ | 2 )ii • x 2 dsids2, (1) 

where Xj(sj) represents the four-position of particle % = 1,2 parametrized by its arc-length 
Si , double bars stand for the four-vector modulus ||a;i — x 2 \\ 2 = {x\ — x 2 ) ■ {x\ — x%), and 
the dot indicates the Minkowski scalar product of four-vectors with the metric tensor g^ v 
(<7oo = 1) 9u — 922 = #33 = — 1) (the speed of light is c = 1). The integration in Eq. ((TJ can 
not be restricted to a segment of trajectory going from an initial to a final time, because 
Lagrangian ((TJ at the end-points involves the future and past light-cones, which are out of the 
segment. This difference from our Poincare- invariant Lagrangian to the usual Hamilton's 
principle of Galilei-invariant mechanics demands that action ([TJ be defined by the whole 
orbit. A sensible mathematical variation of Eq.Q is along trajectory variations 5xi(t) that 
vanish at both t = ±00, such that 5Sp is a sensible finite quantity to be minimized. This 
difference is discussed in and we shall henceforth ignore it and derive our equations 
of motion formally by extremizing action ((TJ. The attractive two-body problem is defined 
by Eq.fTJ with e% = — e 2 = e (Hydrogen atom), while the repulsive problem is defined by 
Eq. (Q) with e\ = e 2 = e. For the electromagnetic two-body problem with arbitrary masses 
Eg. (Ill)'), the only known analytical solution is the circular orbit for the attractive problem 
12L 1 131 . A comprehensive discussion of literature on the Fokker Lagrangian is given in 
Ref. [3j. Some existence and uniqueness results for this two-body system are reviewed in 
Appendix I. 

This paper is divided as follows: In Section II we construct a Hamiltonian that is general 
enough to describe the non-runaway solutions of the Fokker Lagrangian for the two-body 
problem with arbitrary masses and with either attractive or repulsive interaction. In Section 
III we use this Hamiltonian in an EBK quantization procedure that is non-perturbative and 
can provide a set of energies free of infinite quantities. In Section IV we put the discussions 
and conclusion. In Appendix I we construct a Lorentz frame where a generic non-runaway 
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orbit of the Fokker Lagrangian is time-reversible (the CMF frame). 



II. OUTLINE OF THE METHOD 



The method to describe twice-monotonic orbits of the Fokker Lagrangian by an implicit 
Hamiltonian formalism was discussed in Ref. [3J] for the special case of equal masses with 
repulsive interaction. In the following we generalize this result for the case of arbitrary 
masses with either repulsive or attractive interaction. 

n 

The first step of the method consists of a transformation to new variables [11 1 

£ x = ti - xi, Ci = *i + ( 2 ) 
£2 = h -x 2 , C 2 = h + x 2 - 
This transformation splits the action integral into the sum of two separate actions 

S F = l(S a + S b ), (3) 

with 



S a = - J mxid^dQ 1 ' 2 - J m 2 (dt 2 dC 2 ) 1/2 (4) 



-<.<>/ / ?5 X ff KA + ^Ci), 

1^1 — ?2l 



and 



S b = - J mxid^dQ 1 ' 2 - J m 2 {di 2 dC 2 ) 1/2 (5) 
-eiea J J 5 ^ ~ ^ (d^dC 2 + d^dQ. 

Integration of the 5 function in Eq. (@J) gives a nonzero contribution only where £i and 
( 2 take equal values Ci — C2 = C^ an d this ( is the natural evolution parameter of action 
S a (henceforth called type a foliation). Analogously, for Eq. integration over the S 
function produces a nonzero contribution only where £ x = £ 2 = £, and this £ is the natural 
time parameter of action St, (henceforth called type b foliation). Notice that with type a 
foliation particles are in the light-cone condition (xi — x 2 ) 2 — (ti — 1 2 ) 2 = and particle 2 is 
ahead of particle 1 in time by r a . The light-cone distance r a in the type a foliation is 

Ta = -\{^-Q. (6) 



With type b foliation particles are also in the light-cone condition, and particle 2 is behind 
in time by , and the light-cone distance is 

n = ^(Ci-C 2 )- (7) 

After expressing action (j^J) in terms of the time-like parameter (, the associated Euler- 
Lagrange equation is a simple ordinary differential equation. The Euler-Lagrange problem 
for action (JHJ is analogous, with ( replaced by £. A Lagrangian (such as Eqs. (J!J) and ©), 
whose Euler-Lagrange equation is an ordinary differential equation is henceforth called to 
be in the local form. Notice that Lagrangians (0J and (jHJ) are also separately covariant. 

The fact that the Fokker Lagrangian is not in the local form is a motivation to search for 
an equivalent covariant and local Lagrangian, and the natural candidates would be S a and 
St, of Eqs. and (0). With that in mind, one could speculate if an orbit that extremizes 
Eq. (JTJ) could also extremize actions S a (Eq. (J3J)) (5S a = 0) and Sb (Eq. ©) (5S& = 0) 
separately. Each separate minimization, in general, yields a different trajectory, and as 
dlsCUSS ed ln g these are oever ee.ee!. ,n Re, ,3] we found that th e Rectories 

of the following alternative pair of Lagrangians 

SS a = 5G, (8) 

and 

5S b = -5G, (9) 

coincide for a choice of the Lagrangian G. A trajectory that satisfies Eqs. (jHJ) and will 
also extremize the Fokker action (JTJ), a simple consequence of Eqs. (JHJ), (jEJ) and 

5»S F = -5S a + -5S b = -6G - -5G = 0. (10) 
2 2 2 2 

Our problem is then to fix the Lagrangian G such that Eqs. (JBI) and © yield the same 
trajectory and such that Eqs. (jHJ) and Q be ordinary differential equations (rather than 
delay equations). Such G is henceforth called a bilocal ghost Lagrangian. 

In Appendix I we prove that any non-runaway orbit of has a Lorentz frame in 
which it is time- reversible (henceforth called the CMF frame). In the construction that 
follows we operate on this CMF Lorentz frame. Since the orbit is time reversible, we can 
restrict the minimization of Eqs. (jHJ) and Q to the set of time-reversible orbits ([xi(—t) = 
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x\{t)] and [x2(—t) = X2(t) ]). We also add the physical property that both advanced and 
retarded distances decrease monotonically to a point of minimum and then start increasing 
monotonically again. We henceforth call this set of orbits the CMF set. The piecewise 
monotonic property is a consequence of the velocity being a monotonic function of time, 
which was proved in Appendix 1 of Ref. |3| for the repulsive case. The generalization to 
the attractive case along symmetric orbits is trivial. We henceforth refer to a CMF orbit 
simply as a twice- monotonic orbit. Assuming that a twice- monotonic and time- reversible 
orbit exists in our given Lorentz frame, as this orbit extremizes Eq. (0) over all orbits, it is 
necessarily an extremum inside the CMF set. 

In the following we generalize some integral identities first used in Ref. 0] for CMF 
orbits. The time-reversal operation maps type a parametrization onto type b and acts on 
CMF orbits by the following map: £i 2 — > — £ x 2 , £i 2 — > — Ci r a — ► ?"&■ The simplest type 
of integral identity, valid for an arbitrary integrable function <f>{x) of the real variable, is 

4>(r a )d( = / <Kn,R. (ll) 

Jb 

The lower index of the integral denotes the parametrization type, and Eq. ()11|) is a con- 
sequence of the coordinate transformation induced by the time-reversal symmetry of the 
CMF set Q — > — £ (the Jacobian of this transformation is one). For the attractive problem 
we can restrict the integration to twice the finite interval between two consecutive collisions 
(the period), while for the repulsive problem we must define the above integrals over the 
infinite interval. In this work we ignore questions of convergence of the integrals, an am- 
biguity inherited from the Wheeler- Fey nman theory and discussed in Ref. 11]. The same 
time-reversal action (£ 12 — > — £i2i£i2 — > — C1.2) along CMF orbits generates the following 
identities for arbitrary functions Vi(r), V^(r), «i(r) and a^r) of the real variable 

/ai, 2 (»'„)(<if 1 <iCi) 1/2 = /ai,2W(^WCi) 1/2 . (13) 

The above identities suggest that we use a ghost Lagrangian G of type 

1 . . , , • 1 



G = JW a ) + -V x {r a )k x + £ V 2 (r a )£ 2 + ai(r«) + a 2 (r«) V&C- (14) 

The dot over £ x 2 in Eq. (fTlj) indicates derivative respect to C ( the time-parameter of case a) . 
This G is in the local form when added to S a , where £ plays the role of the time parameter 



and the coordinates are £ 1; £ 2 - When this same G is subtracted from action S b , the integral 
identities allow us to express G as 

G = j^{r b ) + ^i(r 6 )C! + \v 2 (r b )C 2 + «i(r 6 )vC + a 2 {r b )sft~ 2 ]d^ (15) 

which is also in the local form for action S b , with £ being the time parameter and the coor- 
dinates being ( x and ( 2 . Notice that the dot over ( 12 in equation (fTKj) indicates derivative 
respect to £ (the time-parameter of case b). Lagrangian (|14j) is the most general ghost 
Lagrangian whose associated Hamiltonian involves quadratic rational functions of the mo- 
menta, and we shall see that it suffices to describe any non-runaway orbit with two monotonic 
branches. 

The condition r = divides the phase space of a twice-monotonic orbit in two disjoint 
regions according to whether r > or r < . In each disjoint region the Lagrangian has a 
branch that can be uniquely inverted to a Hamiltonian formalism. To avoid the overloaded 
notation of Ref. |g, we henceforth drop all indications about branches. Notice that after 
inclusion of the ghost Lagrangian, the problem SS a = SG of Eq. (JHJ) implies the Euler- 
Lagrange equations for L a = S a — G 



L a = - J [ Mla Jl 1 + M 2a ^ 2 + ei _^ 2 + \v^r a % + \v 2 {r a % + <Kr a )]d(, (16) 

where M\ a = mi + «i(r a ) and M 2a = m 2 + a 2 {r a ). The time-reversed problem SS b = —SG 
is described by L b = S b + G, 



L b = -J [M lb ^ 1 + M 2b ^ 2 + | ei _ 62 | - ^i(r 6 )Ci - \v 2 { n )C 2 - <Kr b )]d£, (17) 



ICi-C 

with Mi b = mi ~ <^i{ r b) and M 2b = m 2 — a 2 (r b ) . The Hamiltonian in each case and branch 
is given by 

— 1 /If 2 /If 2 

and 

Notice that the Hamiltonian H a depends only on r a = — §(£i — £ 2 )> sucn that P a = pi +^2 
is a constant of motion. Analogously for type b parametrization, Hamiltonian H b depends 
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only on r& = n(Ci — C2)) sucn that A = Pi + P2 is a constant of motion. The constant 
P a = Pi + P2 of case a suggests the following canonical change of variables 

X = ^(£i + £ 2 ), P = Pi+P2 (20) 

x = ^(£1 -£2), P = Pi -P2- 

For type 6 we use the analogous transformation with £ replaced by £ in Eq. (|2TIj) . We use Eq. 
(j2"UJl to express pi and P2 of Eq. (JT5J) in terms of P and p and substitute into the condition 
H a (r a ,p a , E, P) = E a , yielding a quadratic equation for p, with solutions 



p a {r a} E a ,P a ) = V<2 — + -— =f 



2 + ^ a (£ a + 0) 2 2 r a ; ^(£ a + 0) 2; ' 

(21) 

where Q a = |(^i a + M$ a ) , A a = |(M| a — Mf a ) , r a = |x| and the plus or minus describes 
the two branches. The separation for case b is analogous. Notice that in Ref. [3] we used 
only one function V = \(Vi + V2), a trivial gauge transformation if V\ and V% are functions 
only of r a , but nontrivial if V\ and V2 depend on E and P, as in the generalization that 
follows. This generalization is necessary for the different mass case. 

As we show below, it is always possible to adjust the potentials such that Hamiltonians 
(|T%|) and (fTTjj) have a common trajectory. The drawback is that these so adjusted potentials 
depend on E and P as well. This dependence suggests that we generalize the potentials 
of Eqs. (fTHj) to functions of E and P P|. For example, the potential 4>(r a ) should be 
generalized to = (f>(r a , E a , P a ) in case a and to = (p(r b , E b , P b ) in case b (an analogous 
generalization goes for V\ t 2 and 0:1,2)- The generalized Hamiltonian is still a function of phase 
space, defined implicitly by Eqs. (fT£|) and (fTHj) . because E and P are functions of the phase 
space. This generalization must be done carefully though, because in general the equations 
of motion of Hamiltonians (fTHj) and (fTT?|) with implicitly defined potentials do not correspond 
to the Lagrangian equations of (fTo^) and (JT7j) anymore, but rather to some highly involved 
Lagrangian with a complex dependence on the velocities. On the other hand, simultaneous 
minimization of a pair of Lagrangians with the simple form of (|TH|) and (fTTj) is very desirable, 
because the Wheeler-Feynman equation of motion (SSf = ) is then a consequence of Eq. 
(fTUj) . In the following we construct a pair of auxiliary Lagrangians with the simple form (|TH|l 
and (|17|) for any given orbit. These auxiliary Lagrangians carry different potential functions 
for each orbit. All we need to know about these auxiliary Lagrangians is their existence, 
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as it guarantees that our orbit satisfies the Wheeler- Feynman equation of motion by means 
ofEq. (HUJ). 



The auxiliary Lagrangian is constructed from the potential functions of (r, P, E) as fol- 
lows: For an orbit of energy E Q and momentum P , we define auxiliary fixed form potentials: 
0*(r) = 0(r,P o ,P o ); V£ 2 (r) = Vi >2 (r, E a , P ); a* 2 (r) = a lj2 (r, E , P ), which in turn define 
auxiliary Lagrangians by Eqs. (jlfij) and Q17j) . The orbit of the auxiliary Lagrangian (jlfij) 
satisfies the auxiliary Hamiltonian equation of motion defined by Eq. (JTSJl with the fixed 
form potentials. This equation is in general different from the Hamiltonian equation of 
Eq. (JTSJl with the implicit potentials. To fix this, we construct the potentials such that 
both equations are equal in each case. For this it is sufficient to supplement the following 
conditions in case a 

0, (22) 
0, (23) 



dH a (p,P a ,r, 


E a ) 


dE a 




dH a (p,P a ,r, 


E a ) 


dP a 




dH b (p,P b ,r, 


E b ) 


dE b 




dH b (p } P b} r } 


E a ) 


dP 



and analogously in case b 

0, (24) 
0. (25) 

These are four conditions imposed on the five potentials 0, V\p and a\^- Before we proceed 
to construct potentials satisfying the above conditions, it is convenient to find the explicit 
solution of Hamiltonian (|18|). which is accomplished by use of a canonical transformation 
with a generating function S given by 

S = PX + W(x, P, E) - E( } (26) 

where W(x, P, E) is defined to satisfy the condition p = dW/dx , with p given by Eq. (|2*T|) . 
The new momentum associated with the old variable X is the same old constant P = dS/ dX 
and the other new momentum is the energy E (with this last definition we exploit the fact 
that E is already one argument of the implicit potentials). We choose S in the manner 
of Hamilton- Jacobi such that the new Hamiltonian vanishes: K = H + = 0. As the 
Hamiltonian is zero, the new coordinates are defined simply by two constants X and Cq 
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X = dS/dP = X + dW/dP (27) 
C = -dS/dE = C - dW/dE 

Eq. (|2~7j) defines ( and X as functions of r a = |x| for case a (and £ and X as a function of 
ri, = \x\ for case b ). We shall also need the differential form of Eq. (|27|) relative to x, 

dX = -(d 2 W/dxdP)dx = -(dp/dP)dx, (28) 
d( = (dW/dxdE)dx = (dp/dE)dx } (29) 

where we have used p = dW/dx (definition of the Hamilton- Jacobi transformation) and 
exchanged the partial derivatives. The differentials dxi and dti along the trajectory are 
obtained using Eq. (J2J) to relate particle coordinates to X and Q followed by use of Eqs. 
(|2*H|) and (j2HJ) to relate dX and d( to dr a = \dx\ . The explicit solution is 



dtla,b = 


2 [ dP 


dp a ,b 1W 

dE l)dT ^ 


(30) 


dt2a,b = 


2 [ dp 


QE +l)dr a>b , 




dXi a ^b = 


2 [ dP 


dp a ,b N , 

dE 




dX2a,b = 


2 [ dP 


^ +l)rfr a , b , 





where the plus sign holds for case a and the minus holds for case b in the time variables. 
Eq. (plj) gives the explicit solution in terms of p a ,b( r , E, P), as given by Eq. (|2*Tj). 

As time reversal maps case a into case b, a time-reversible orbit of both (jl8J) and (fT9"|) 
must satisfy v\ a {r) = —Vib(r) and V2 a (r) = —v 2b {r) for r e [r c ,oo]. An economical way to 
satisfy this condition is to express the functions dp/dP and dp/dE as 



dp a dp b _ cosh[.s(r,P,P)] 

dP ~ dP ~ sinh[s(r,P)] ' 1 j 

dPa = dPb= F(r,E,P) 

dE " «9P ~ sinh[s(r, P, P)] ' 1 ' 

such that Ui a (r) = — t>i&(r) and t>2a(?") = — i>2b(r) are automatically satisfied by use of Eqs. 
(jnilj). We are now ready to obtain an analytical solution of Eqs. (|2^j) . (J2I|), (|3*Tj) and (|3^j) in 
terms of the undetermined functions s(r, P, P) and F(s, E, P) defined by Eqs. (}3"Tj) and 
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Our solution is defined such that E a = Eb = E and P a = Pb = P- For case a the equations 
of motion for £ x and £ 2 derived from Hamiltonian ( ITHj) with use of conditions are 

exp(s) (mi + ol\) 2 



and 



d( ' F(r) ~ 4(p 1 + §F 1 + ^)2 : 
d£ 2 exp(-s) (m 2 + a 2 ) 2 



(33) 



(34) 



rfC F(r) ~ A{p 2 + \V 2 + e ^y 
As discussed in Ref.j^J, there are four possible square roots of Eqs. (jSHj) an d (EH), but 
only one choice does not pose a constraint involving s(r, E, P) and F(r, E, P), leaving both 
functions arbitrary. This physical choice of signs is defined by 

P + Va + V 1 + — = -VF(m l + a 1 )exp(-s/2) 7 (35) 
r 

P-p a + V 2 + — = v / F(m 2 + a 2 )exp(s/2). (36) 

where we have used p\ = \{P a + Pa) and p 2 = \{P a — Pa) with P a = P. Case b is analogous, 
and the motion of £i and ( 2 derived from Hamiltonian ( IT9")) with condition Eq. (|24|) defines 
the physical choice of the square roots by 



P + p h -Vi + — = V^(mi-ai)exp(-s/2), (37) 
eie 2 



ei£2 
r 

P-p b -V 2 + ^ = -v / F(m 2 -« 2 )exp( S /2). (38) 
We can use Eqs. (I35j) and (}36|) to eliminate pi and p 2 from Eq. (jTHj) with E a = E, yielding 

E + <p = ——[(rm + a x ) exp(s/2) - (m 2 + a 2 ) exp(-s/2)], (39) 
2y F 

and an analogous use of Eqs. (pTTj). (J35j) and (fT9"|) with E b = E yields 

E-<p = W(mi - ai) exp(s/2) - (m 2 - a 2 ) exp(-s/2)]. (40) 

2v .F 



Notice that Eqs. (j3Tj) and ([32)1 are immediately satisfied if we assume 

Pa = Pb = p(E,P,r). (41) 

Eqs. (jBTj) and (J32*j) actually allow p a and p& to differ by a trivial Gauge function g(r) entering 
in Eq. (|41|) . but that can be absorbed in the definition of Vi and V 2 . We can satisfy the seven 



11 



conditions of Eqs. ()35|) - ()41|) with completely arbitrary functions F(E,P,r) and s(E,P,r) 
in the following solution 

PPexp( S /2) + (P + a^)exp(-s/2) , 

V r sinh(s) 
PPexp(-s/2" 

«2 



PPexp(-s/2) + (P + sia) exp(s/2) 



Psinh(s) 
Vi = —mivF exp(— s/2), 

V 2 = m 2 v / Fexp(s/2), 

mx exp(s/2) — m 2 exp(— s/2) 



and last with p given by 



2^ 

[EF + (P + ) cosh( g )] 

P= r-r-^r . 43 

sinn(s) 

Now that we have solved simultaneously all conditions (jSU), (J32J) and (fSKj) - ^ . 

some remarks are in order: 

(i) To see the generality of our method, we can start from a numerical solution of the 
delay equation of motion, with energy E, and simply adjust the two arbitrary functions 
s(E, P, r) and F(E, P, r) such that Hamiltonian (fT%j) has this orbit as solution. This can be 
done with Eqs. and ()34j) and takes exactly two functions, to describe the two particle 
positions as a function of the light-cone separation. The method has a branch point at s — 0, 
such that for each separation there should be two values for the position, one in the inbound 
flight and another in the outbound flight, such that we can describe twice- mo notonic orbits 
only. A general proof of the twice-monotonic property is given in Appendix A of Ref. 0. 
At the moment the only available method to solve the delay equations of motion is by 
numerical calculation. An automatic consequence of (pT7jl and (J3%|) is that the orbit also 
satisfies Hamitonian ()19p. 

(ii) Condition was used as an ingredient in obtaining Eqs. (l3*3*jl and (|3*4*jl of case a, 
but we never forced it directly. An instructive direct check is to take the partial derivative 
respect to E of the implicit energy equation (E = H a (r,p, E, P) ) defined by Eq. fl~8|) . 
together with substitution of the potentials of Eqs. (jj2j) and use of Eq. (JH2J)- The result is 
dH a /dE = 0, and analogously for case b we obtain 8H a /dE = 0, which are the upper lines 
of (j22j) an d (|24|). The lower lines of (}2*2*|) and ()24j) are obtained in the same way, by taking 
the derivatives of the implicit energy equation respect to P in each case. 
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(iii) After we determine s(E,P,r) and F(E,P,r) , substitution of E = E and P = P Q 
into the potentials produces auxiliary Lagrangians that have an orbit in common, precisely 
because of (ii). And of course, if we did not know that yet, we would prove that this orbit 
satisfies the Wheeler- Fey nman equation of motion (SSp = 0), as a consequence of Eq. (jlUj) . 

Last, as discussed in Ref. [3J, at the branch point s = and r = r Q the continuity 
condition for the functions of Eqs. ()42|) and (|43|) calculates P as 

P = —EF(r ) - (44) 

which will be used in the next section. 

III. INFINITE-FREE EBK QUANTIZATION AND COVARIANCE 



Under a Lorentz transformation with boost parameter w the light-cone distance of case 
a transforms like 



r a = sj{l + w)/{l-w)r a , (45) 

while in case b the same formula holds with the replacement w —>■ —w. Equation (|43j) shows 
that the twice-monotonic property is a Lorentz-invariant property, suggesting that there is 
no obstacle to describe a CMF orbit in another Lorentz frame. 

The Lorentz transformation acts on the coordinates of Hamiltonian (fTH|) (the £ x and £ 2 
defined in Eq. ( |21 ) ) as a simple rescaling 

li = Ui, (46) 

£2 = T - £2 

and also rescales the evolution parameter ( of case a to 

C = A a C, (47) 

with A a given by 
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Substituting the above Lorentz-transformed quantities into (|16|) we find another Lagrangian 
with the same functional form, which proves that a CMF orbit is described in any Lorentz 
frame by a Lagrangian L a of the form To see how the potentials are transformed in 

the new frame, we make use of the explicit solution (J3U|) in both frames, together with the 
condition that the new coordinates Xi a ,ii a , xu and in, must be the Lorentz transformed of 
the old coordinates xi a ,ti a , x\b and ti^with the boost parameter w, yielding 

dp a (1 - w) dp a 



<)i-: ~ (i • tr)oi: 

dpg dp a 

dP ~ dP' 



(49) 



where we used the differential version of Eq. (|43j) ( dr a = J (1 + w)/(l — w)dr a ). Equation 
fl4T?|) shows that s a (r) as defined by Eqs. (pIT]) and (J3~2~)) is a Lorentz scalar: 

s a (f a ,E,P) = s a (X a f a ,E,P), (50) 

and that F a {r a ) (as defined by Eqs. (j3"Tj) and (J3"2~j) ) transforms like 

F a (f a ,E,P) = X 2 a F(X a f a ,E,P). (51) 

Case b transforms in the same way with = 1/X a , which is Eq. f!48|) with the substitution 
w — > —w. This last equation allows us to express the Hamiltonian in any frame by use of 
the CMF form of F(r) and a rescaling depending on the boost parameter. 

It can be seen that the Lorentz transformation of Eqs. (|4*E|) is a canonical transformation 
that scales the momenta with the inverse factor, pi = X a p\ and p2 = X a p2, such that the 
canonically transformed Hamiltonian is 

4 (Pi + l 2 V i + (ft + 1^2 + ^) 

with 0(f a ) = ±<f>(r a ), V(f a ) = X a V(r a ) , Ml(f a ) = M* a (r a ), and M| a (f a ) = M| a (r a ). 
Notice that the new Hamiltonian H a picked a multiplicative factor of A a , and if we perform a 
change to the natural evolution parameter ( = X a ( of the new Lorentz frame, it compensates 
exactly for that factor, going back to the form ()18j) . such that the Hamiltonian has the same 
form in all Lorentz frames (only the value of the energy is changed to E = E/X a ). 

As the Lorentz transformation is canonical for our I hamiltonian, the action / = <j> pdr is 
a Lorentz- invariant, such that EBK quantization 



23, 



24j is a sensible covariant quantization 
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procedure. The need to follow the canonical transformation by a rescaling can be avoided 
by use of a Lorentz- invariant evolution parameter in the Lagrangian of Eq. (£Q). For example 
a Lorentz- invariant parameter can be constructed with the help of Eqs. (J51|) and (J47j) as 

du a = -, ^ = , ^ (53) 
^F a (r a ) y/F a (f a ) 

A 8 the _ , . M/at - .dependent of pa—ion, the action h,e graI 
marks the same orbits with or without rescaling. 

The action Im f+ for the active prooie. <„« can b e „d, h y use of 

Eqs. (g5D . (jU D andl32l as 

/ = -/[^ + ^y(P-^)K, (54) 

and substitution of P as given by Eq. (jUj) into Eq. (J53J) yields 

1 F(V ) 

m - (55) 



rF(r) F(r) 

The function F(r a ) can be calculated by elimination of s(r) from Eqs. fpTTTf) and and 
use of (J2J) as 

^(r a ) = a + + (56) 

where v% is the velocity of particle 1 and t>2 a is the velocity of particle 2 in the advanced 
light-cone parametrized by r a . For a shallow energy orbit, the particles have a low velocity 
in a large portion of the orbit, such that Eq. (|56[) predicts F(r) ~ 1. The main difference 
from Coulombian orbits to our covariant orbits at shallow energies is at the collision, where 
Coulombian orbits reach an infinite velocity, while the relativistic motion saturates at the 
speed of light. Disregarding the collision region, substitution of F(r) = 1 together with the 
low- velocity approximation d( ~ dt (valid away from the collision), into Eq. ([55)1 . yields 
exactly the one-dimensional Coulombian WKB integral 

Equation ()57|) is easily derived using the Coulombian energy and the Levi-Civita regular- 
ization, see for example the Appendix on Levi-Civita regularization of Ref. Q|. At this 
point the reader can see that our strange looking Hamiltonian formalism actually reduces 
nicely to the Coulomb approximation for shallow orbits. The use of / = {n + together 
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with the approximate formula ()57j) predicts the well known Balmer energy terms for the 
one-dimensional Coulombian two-body system. The corrections can be evaluated using the 



exact covariant formula (1551) along the numerical orbits, as obtained in Ref. 19] for the 

n n 

equal-mass case. As discussed in Ref. [19] for the equal-mass problem |1 91] . at the outbound 
collision F(r a ) is singular like l/y/r, and d( = (1 + v\)dt ~ v2r a dr a , which makes the inte- 
gral of (|5*5*j) finite (actually, for the equal-mass case the collision does not even contribute to 
the integral). The merit of this procedure is that it calculates a Lamb shift without infinities, 
which might prove useful for understanding the renormalization of QED. For example, on 
the equal-mass two-body problem, the orbits can be calculated numerically using the regular 
delay equation developed in Ref. . Along those regular collision orbits, there is nothing 
infinite involved in the evaluation of the action integral of Eq. 1)5 5 j) and the above covariant 
EBK quantization procedure yields an electrodynamics free of infinities for the positron- 
ium atom (the quantum energies are all finite). It would be highly desirable to generalize 
this procedure to the 3-dimensional motion of the electromagnetic two-body problem, but 
at present we do not know how to do that. The EBK quantization of Wheeler- Feynman 
theory is discussed in j^j, |3| along with the only physical result so far (i.e. simple EBK 
quantization of circular orbits). This quantization predicts energies with finite small shifts 
in agreement with the Dirac equation to order a 4 j^. This is exactly what one expects 
physically, as those are circular orbits with a finite angular momentum. The infinite Lamb 
shift plagues only zero angular-momentum states, which corresponds to the orbits studied 
in the present work. The possibility of using the regularization of Ref. to produce an 
electrodynamics free of infinities for the positronium atom suggests that further efforts to 
understand the regularization of the different-mass electromagnetic two-body problem could 
shine light on the renormalization of QED. 

IV. CONCLUSIONS AND DISCUSSION 



The idea to remove the field degrees of freedom from Maxwell's electrodynamics goes 
back to Dirac |8| and later Wheeler and Feynman planned to quantize WF2B as a means to 
avoid the divergencies of QED, as in the action-at-a-distance theory the infinite number of 
field degrees of freedom are absent. The difficulties in converting the Fokker Lagrangian 
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;o Hamiltonian form have caused the famous seminar that never came from Wheeler (see 
lo| . page 97). In chapter 5, page 97 of reference Feynman says that ' I didn 't solve it 
either — a quantum theory of half-advanced half-retarded potentials — and I worked on it for 
years... '. This is still an outstanding problem today and the difficulties in casting relativistic 



28, 



Lagrangian interactions into Hamiltonian form are explained in references 

Our description does not violate the no- interaction theorem , and there are two places 
where we avoid it: (i) the no-interaction theorem is an obstacle to covariant Hamiltonian 
description of two interacting particles only in fully 3-dimensional motion. We are restricted 
to 1-dimensional motion; (ii) Because the evolution parameter in our Hamiltonians is not 
time but rather ( for case a (and £ for case b), the no-interaction theorem does not apply. 
In principle, because time is not the evolution parameter, even in three dimensions the no- 
interaction theorem would not be an obstacle to a Hamiltonian description, and that is an 
open problem. The merit of our constructive approach is that it starts from the physical 
theory of action-at-a-distance and does not involve expansions. This procedure should be 
compared to the very interesting covariant extrapolation of Wheeler- Feynman theory within 
constraint dynamics, which is a covariant way to introduce the correction of order {v/c) 2 to 
the Coulomb problem j^. 

The energy equation for our Hamiltonian ()18|) is the following quadratic equation for p : 

p 2 + A(E,r)p + U(E,r) = 0, (58) 



with 



A(E,r) = [ ^7^ + (V 2 ~ Vi)), (59) 

jj( F r ) = Ko+^L) /p | eiea (m 2 a V 2 + m 2 a V x ) 
1 ' ] 2(E + <f)) 1 + r ' 2(E + <P) " { } 

where m la = mi + a%(r) , m 2a = rr^ + a^?") , V\(r) and Vj>(r) are determined from s(r, P, E) 
and F(r,P,E) by Eq. The procedure outlined above determines the potentials only 

along the orbit, but it would be desirable to extend these potentials to the whole phase 
with some analytical continuation and/or use of symmetry. If that is accomplished, after a 
symmetrization, Eq. (|H%|) can be used in a quantum canonical formalism with p and P being 
the quantum momentum operators. The trivial separation of the center-of-mass coordinate 
X produces a familiar looking Schroedinger equation. This equation should be compared 

17 



and 



to the one-dimensional Klein-Gordon equation, that has been studied in connection with 

n 

the one-dimensional hydrogen atom [25||. The quantization of the one-dimensional hydrogen 
atom interacting with the electromagnetic fields is an interesting testbed problem, as the 
minimal coupling to the electromagnetic field produces divergencies in perturbation theory 
that are exactly like in the 3-dimensional problem (and which are absent in the quantization 
outlined above). In that respect, Eq. (J55j) is a properly renormalized version of the Coulomb 
Hamiltonian operator. 

Last, our procedure could be applied to the problem of calculating level shifts for an atom 
in the presence of a strong gravitational field, by incorporating the metric tensor g^ v in the 



definition of the variables of Eq. (J2J), which should be modified to ^ l = y/g^ti — J\gn\xi 
and Ci = Jgooh + Jgrixi. 



V. APPENDIX I: CONSTRUCTION OF THE CMF FRAME 



For the repulsive two-electron problem (mi = 1712 = m and e\ = €2 = e) and along 
symmetric orbits [— ^(i) = x i{t) = x (t)], the minimization of action (Q) yields the following 
equation of motion 

d v e 2 / l-v(t-r) \ e 2 / 1 + v(t + q) \ 

m Jt { ^T^2 } - 2^ \i + v (t-r) ) + V [i- v (t + q) ) ' { } 

where v(t) = dx/dt is the velocity of the first electron, of mass m and charge e, and the 
functions r(t) and q(t) are the time-dependent delay and advance, defined implicitly by the 
light-cone conditions 

r(t) = x(t)+x(t-r), (62) 
q(t) = x{t)+x(t + q). 

In general, a neutral-delay equation such as Eq. (JBTj) requires a pair of world-line segments 
as the initial condition (one world-line segment for each particle). As discussed in Ref. 

, the initial world-line segments can be provided in such a way that Eq.([61|) is well- 
posed by using "maximal independent segments ". A pair of world-line segments is called 
independent if the end points of each segment lie on the forward and backward light-cones 
of a single point interior to the other segment. A surprising existence theorem was proved 
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in Ref.([15j) for the problem of Eq. (j61|) along symmetric orbits [—X2(t) = x±(t) = x(t) ] 
using a Banach to Banach contraction mapping; This theorem states that for sufficiently 
low energies, Newtonian initial conditions [ x(0) = x Q and v(0) = v ] determine the unique 
symmetric orbit that is globally defined (i.e., that does not runaway at some point) 15 1. 

Later in 1981 it was proved for this same equal- mass problem that the delay Eq. (J61I ) is 
equivalent to a local Newtonian equation |lfij |. 

The physical meaning of the above order-reduction of the delay equation is the following: 
of the infinite many initial conditions necessary to solve the delay equation of motion, only 
the Newtonian type "position and velocity values at t — " suffice to determine the non- 
runaway orbit. The remaining conditions either take values determined by the initial position 
and velocity or the orbit will be a run-away. The above existence/uniqueness results are 
certainly true in the low-velocity Coulombian limit, as the Coulomb ODE's have unique 
solutions determined by initial position and velocity. Therefore the existence/uniqueness 
result is either true in general or at least up to some high energy where strange new solutions 
bifurcate out. In the following we use this Newtonian order- reduction to define the CMF 
frame. 

We consider here the two-body problem of Eq. (|TJ with attractive interaction and arbi- 
trary masses. The construction of the Hamiltonian given in Section II for time-reversible 
orbits can be generalized to arbitrary orbits if we are able to find a Lorentz frame where 
this generic non-runaway orbit is time-reversible. This frame is henceforth called the CMF 
frame for the orbit. To accomplish this we start below by proving the following Lemma: 

Lemma 1: Given any non-runaway orbit of Eq. (JljL defined in a given Lorentz frame, we 
can construct another Lorentz frame where the two particles come to rest at the same time 
(t = 0) along this orbit. 

The proof will be given for the attractive case and we define the measure time such that 
at t — the particles are at the origin and moving in opposite directions with the speed of 
light: xi(0) = £2(0) = 0, xi(0) = 1 and ±2(0) = — 1 (the outbound collision). As the force is 
always attractive, the absolute velocity of each particle will decrease monotonically until it 
comes to rest at some turning time, and then it accelerates again to the speed of light until 
the next collision, that happens at a later time t = T (the inbound collision). At this second 
collision the particles bounce at each other and start again with outbound velocities. The 
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proof can be made as well for the case where particles pass through each other at the inbound 
collision, such that the orbit is twice as long. For simplicity we assume that they bounce, 
such that the unit cell of our periodic orbit is the interval (0,T), and particle one is always 
on the right-hand side. The turning times of the particles are in general two different values 
in the interval (0, T). To fix that, we use a Lorentz transformation with boost parameter 
w, defined such that the transformed of events (ti,xx(ti)) and (t 2 ,x 2 (t 2 )) be simultaneous 
in the new frame. This imposes the condition 

ti — wxi(ti)=t2—wx 2 (t 2 ). (63) 

As we also want that both particles be at rest simultaneously in the new frame, we must 
choose w such that 

xi(ti) = x 2 {t 2 ) = w. (64) 

Because the velocities are monotonic in the interval (0, T), condition (jHlj) defines t 2 as a one- 
to-one function t 2 (t\) from the interval (0,T) onto itself. We can then rewrite condition 
(|63|) as the zero of a function of the variable ti of the interval (0, T) 

F(h) =h- t 2 (h) - - x 2 (t 2 )) = 0. (65) 

Using the facts that t 2 (ti = 0) = T, ii(0) = land t 2 (ti = T) = 0, X\(T) = — 1, and also that 
^1(^1) — £2(^2) > for all times, we find that the continuous function F(t{) is negative at 
t\ = and positive at ti = T, such that it must have a root in the interval (0, T). This value 
of t\ calculates the boost w by Eq. (jfi4j) and defines a Lorentz frame where both particles 
are at rest simultaneously (the CMF), as we wanted to demonstrate. 

We have then constructed the CMF frame with the property that both particles have 
zero velocity at the same time, say t = 0, which we can now use as "Newtonian initial 
conditions " . As discussed above, position and velocity initial conditions determine a unique 
non- runaway orbit, and since an initial condition with both particles simultaneously at rest 
is invariant by time reversal, the unique orbit determined by it must be time-reversible. 
We actually need the generalization of the Banach to Banach proof of Refs. (^^) for the 
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attractive case with different masses to complete our proof. 
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